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I. INTRODUCTION 

The influence that cold front passages have on Louisiana coastal environments, 
including land loss and land building processes, has been the primary topic of this 
multidisciplinary research. This research has combined meteorological, remote sensing, 
and coastal expertise from the University of Wisconsin (UW) and Louisiana State 
University (LSU). Analyzed data sets include remotely sensed radiometric data (AVHRR 
on NOAA-12,13,14, Multispectral Atmospheric Mapping Sensor (MAMS) and MODIS 
Airborne Simulator (MAS) on NASA ER-2), U.S. Army Corps of Engineers (USACE) 
water level data, water quality data from the Coastal Studies Institute (CSI) at LSU, 
USACE river discharge data, National Weather Service (NWS) and CSI wind in situ 
measurements, geomorphic measurements from aerial photography (NASA ER-2 and 
Learjet), and CSI ground based sediment burial pipes (for monitoring topographic change 
along the Louisiana coast) and sediment cores. The work reported here-in is a 
continuation of an initial investigation into coastal Louisiana landform modification by 
cold front systems. That initial effort demonstrated the importance of cold front winds in 
the Atchafalaya Bay sediment plume distribution (Moeller et al. 1993), documented the 
sediment transport and deposition process of the western Louisiana coast (Huh et al. 

1991) and developed tools (e g. water types identification, suspended solids estimation) 
from multispectral radiometric data for application to the current study. This study has 
extended that work, developing a Geomorphic Impact Index (GI 2 ) for relating 
atmospheric forcing to coastal response and new tools to measure water motion and 
sediment transport. 

Cold fronts impact coastal regions through their associated wind, temperature, and 
pressure patterns, causing redistribution of suspended sediment (SSC), water level 
fluctuation, wave action, and by affecting water temperatures and moisture content of 
fresh shoreline sediment deposits (Roberts et al. 1987). This study has contributed 
analyses of coastal circulation, suspended sediment transport, water types, and 
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geomorphic response to wind system forcing by cold front passages. Primary findings of 

this study are: 

1. Geomorphic impact on Louisiana’s muddy coast region by cold front passages 
can be parameterized using wind energy in combination with water level. In 
particular, only during high water level events is the cold front wind energy 
effective in imparting geomorphic impact. This finding is the basis of the 
Geomorphic Impact Index (GI 2 ), which can be used to assess seasonal to 
interannual impact on coastal landforms. 

2. Estimates of water motion and sediment transport indicate large short term 
variation of nearshore circulation and mass transport, intermittently feeding and 
starving coastal estuaries and nearshore waters with land-building material. 

The implications of this finding for maintaining coastal marine ecosystems are 
not yet fully understood. This finding was made possible by the use of the NASA 
ER-2 research aircraft, which had the capability to collect high spatial (100m or 
better) and temporal (1-2 hours) resolution data. 

3. A method of multispectral coastal water type analysis has been developed using a 
combination of water surface temperature and visible/near IR subsurface 
reflectance in a postfrontal scatter plot analysis. Remote sensing measurements 
made by MAMS (100m resolution) and MAS (50m resolution) on the ER-2 show 
five distinguishable water types including river water, estuarine water, inner 
shelf water, marsh drainage, and soil water drainage (Huh et al. 1996) in the 
Atchafalaya Bay area of Louisiana. The water types illustrate coastal processes 
such as marsh drainage, soil de-watering, and estuarine-saltwater exchange. 

4. Frequent, synoptic scale (but lower energy) cold front systems appear to be more 
important than the infrequent, generally localized (but powerful) tropical 
systems. However, longer periods of record are necessary to fully assess impacts 
of cold fronts and tropical systems. 
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Because of the difficulty of gathering surface data in inaccessible coastal 
environments, this study was designed to apply remote sensing as a primary source of 
information for coastal research. This project thus served as a testing ground for the 
informational content of high resolution remote sensing measurements in a coastal 
environment. The understandings gained from this effort will benefit the application of 
multispectral observations from the MODIS instrument (estimated launch in 1998). 

n. RESEARCH ACCOMPLISHMENTS 

A. Data Collection and Observations 

1 . Remote Sensing and Ancillary Data Collection 

MAMS and, starting in 1 994, MAS were flown on NASA's ER-2 aircraft over the 
Louisiana coast beginning with the cold front season of 1987-88. A detailed description 
of the MAMS instrument is available in Jedlovec et al. (1989); details on MAS can be 
found in King et al. (1996). A list of Louisiana coast flights through the end of this grant 
is shown in Table 1. The ER-2 flights comprise an eight year, high spatial resolution 
remote sensing data record of the Louisiana coastal region. The data flights were planned 
to optimize clear sky data collection over the coastal zone. During most flights, CSI boat 
teams collected water quality measurements in the Atchafalaya Bay region. In situ 
measurements included suspended sediment concentration (SSC), sea surface temperature 
(SST), secchi disk depth (SDD), and CTD measurements. The MAMS and MAS data 
sets have been split into straight line flight track segments, navigated, and archived on tape 
for future use. 

In addition, the NASA ER-2 aircraft collected 5 and 10 meter resolution false 
color infrared aerial photography data during most flights. These data were used to 
monitor small scale changes in geomorphology, and also served as a verification tool for 
the MAMS and MAS observations. 
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AVHRR imagery from the NOAA series polar orbiting satellites was used to 
maintain a daily surveillance of the turbid waters of the Atchafalaya Bay/Chenier Plain 
coastal sedimentary system. Both sea surface temperature and suspended sediment 
calculations were made from the AVHRR 1 km imagery at the LSU Earth Scan 


Table 1. - MAMS/MAS Louisiana Coast Flight Summary 


DATE 1 

COVERAGE 2 

in situ 

COMMENTS 

01/27/88 

Type 1 


postfrontal, low water 

08/26/88 

Type 1 


summer; thin cirrus 

11/02/88 

Type 1 

X 

cirrus 

03/30/89 

Type 1 


prefrontal; SW wind 




high water 

04/01/89 

Type 1 

X 

postfrontal; NE wind; 




FROPA pair w/Mar30 

04/06/89 

Type 2 


triple overpass; high water 




transition to prefrontal 

12/13/89 

Type 1 

At 

postfrontal; low water 

04/15/90 

Type 2 


triple overpass; postfrontal; 




high water 

04/16/90 

Type 1 


transition to prefrontal; high 




water 

12/04/90 

Type 2 

X 

postfrontal; low water 

12/05/90 

Type 1 

X 

postfrontal; light winds 




consecutive w/Dec 4 

03/30/91 

Type 2 


triple overpass; postfrontal; 




cloud streets 

04/01/91 

Type 1 


postfrontal 

11/24/91 

Type 1 


MAS flyby 

12/11/92 

Type 1 

X 

postfrontal; some cirrus 

02/02/94 

Type 1 

X 

MAS; transition to prefrontal 

02/14/94 

Type 2 

X 

MAS; triple pass; postfrontal 

02/17/94 

Type 1 

X 

MAS; postfrontal; stratus; 




high water 

01/08/95 

Type 1 

X 

MAS; postfrontal; high water 

01/24/95 

Type 2 

X 

MAS; postfrontal; triple pass 

04/09/96 

Type 2 

X 

MAS; postfrontal; triple pass 




high water 

Notes: 


1 For flights before 1995, 12 channel digitizer used (7 channels @ 8 bit precision, 

4 channels @ 10 bit precision, 1 channel for extra bit bucket). Beginning in Jan 
1995, 50 channel digitizer used. (50 channels @ 12 bit) 

2. Type 1 flights cover full Louisiana coastline in 4 flight lines. Type 2 flights use 
racetrack pattern to repeatedly overfly selected sites (typically the Atchafalaya Bay 
region of the Louisiana coast). 

* Boat ramp water level too low to launch boat. 
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Laboratory, a multisatellite ground station. A suspended sediment algorithm (Stumpf, 

1 992) using Rayleigh corrected channel 1 minus channel 2 was used to map the 
distribution of turbid waters. The algorithm has been calibrated against water samples 
collected by helicopter (Walker et al. 1996). This allowed detection of the size and 
orientation of river/estuarine discharges, and flow paths of the suspended sediments 
throughout the region. The AVHRR imagery was also used to direct field activities for 
surface in situ sampling. 

In addition to remote sensing and boat in situ data sets, the following data sets 
were collected: 

1) National Weather Service surface wind data for the years 1990-1997 at Lake Charles 
and New Orleans, LA. 

2) Atchafalaya River discharge (daily) and sediment load (bi-weekly) for 1980-1997. 

3) Water level data at the Amerada Hess platform near Eugene Island, 1987-1997. 

4) Special surface wind (1995-1997) and tidal gauge (July 1995) data sets at instrumented 
Fina #13 platform in nearshore Gulf of Mexico for 1995-1997. 

5) Deployment of topographic sediment monitoring pipes and seasonal measurements of 
topographic change (1994-1997). 

These data were applied primarily to the formulation and analysis of the Geomorphic 
Impact Index, which relates meteorological forcing of cold fronts to coastal 
geomorphology. 

2. Suspended Sediment Monitoring and Transport 

Suspended sediment distribution and transport were monitored using MAMS and 
MAS atmospherically corrected visible/near IR reflectances, calibrated to SSC using in 
situ measurements (from boat teams deployed by CSI during ER-2 flights). Atmospheric 
correction was based on work by Gordon et al. (1983), Guzzi and Zibordi (1987), and 
Gumley et al. (1990). Repeat coverage MAS and MAMS SSC maps proved useful for 
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estimating water motion and sediment transport in the turbid Atchafalaya Bay region 
(Moeller et al. 1996). Maps of water velocity and sediment transport demonstrate 
significant small scale variability and show that velocity and sediment transport are not 
highly correlated (Figure 1). Transport of suspended material is central to both terrestrial 
and subaqueous geomorphic evolution. Effects of sediment deposition on the delicate 
estuarine ecosystem are not well understood at this time but are of great economic interest 
(e g. destruction of oyster beds, freshwater versus brackish water domination). 

3. Coastal Water Types Detection 

Coastal waters in plan view are generally a mosaic of waters of different 
temperatures and turbidities (Huh et al. 1996). Water types are particularly variant in 
regions such as the Louisiana coast where the Mississippi River and estuarine discharges 
create a range of suspended sediment concentrations ( 0-1200mg/l). Experimentation with 
MAMS imagery resulted in the critical observation that there is an inverse correlation 
between the visible/near infrared reflectances and the water surface temperature (Moeller 
et al. 1989). This correlation made possible a temperature/turbidity analysis similar to the 
classical temperature/salinity analysis essential to identifying oceanic water types. Water 
type analysis has turned out to be a useful tool in the study of coastal circulation and water 
type formation processes, particularly after a cold front passage. 

B. Geomorphic Developments in the S.W. Louisiana Coastal Zone 

Investigation into relationships between external influences (e.g. cold front wind 
systems, water level variation, etc.) and geomorphology of the Louisiana coast were a 
primary topic of the project. Specifically, it has been observed that the easternmost portion 
of the Louisiana Chenier Plain has changed within the last decade from an erosional to a 
progradational sedimentary regime. The progradational/accretional deposits are sheets of 
fine grained sediment. Geological studies of these washover deposits implicate the 


6 




Figure 1 . MAS near infrared imagery of Atchafalaya Bay region of the central Louisiana 
coast on 24 January 1995. Image is remapped from aircraft projection to Mercator 
projection. On 24 January, the ER-2 Hew three consecutive overpasses of the Atchafalaya 
Bay in a two hour period. The images were remapped to Mercator projection, animated, 
and targets were selected and tracked from image to image providing the water motion 
(green arrows) analysis shown. Sediment transport estimates (red arrows) were produced 
by combining water motion vectors with MAS SSC map. Water motion errors due to 
“jitter” in remapped imagery are estimated at about 10-20# of vector speed (based on 
tracking land jitter). MAS SSC regression uncertainty of 38mg/l is also about 10-20# 
error. 






southerly winds and storm surges of cold fronts and other storms as primary forcing 
mechanisms for coastal accretion and progradation. The approach chosen was to measure 
the changes in the coastal deposits at the particularly active sites of western Louisiana and 
relate those measurements to the external forcing. The ultimate goal of this effort was to 
develop a Geomorphic Impact Index (GI 2 ), which amalgamates the external influences into 
a measure of actual or potential coastal landform change. 

1 . Geomorphological Surveys and Measurements 

Measurements of lateral accretion, vertical accretion, and erosion along the 
western Louisiana coast were made using aerial photography (NASA ER-2 and Learjet) 
and arrays of ground based sediment burial pipes. Aerial photography (1987-1997) was 
used to estimate lateral changes both shore normal and shore parallel; an array of sediment 
burial pipes (1994-1997) were used to monitor the vertical accretion. The sediment burial 
pipes remain in place as useful reference points from which to monitor geomorphic change 
in years to come. 

Land accretion and erosion in the aerial photography was measured manually using 
identifiable landmarks in the photography at several sites along the western Louisiana 
coast. Of primary interest is the muddy prograding coastal region westward (and 
downdrift) of the sediment rich Atchafalaya River mouth and Bay. Identifiable landmarks 
(canals, structures, etc.) along this muddy coast were used as “anchor” sites to measure 
mud area over a 10 year period of aerial photographic surveys (some 40 different dates, 
from 1987-1997). Active erosional sites to the west and east of the muddy coast were 
also monitored. 

Sediment burial pipes (galvanized iron) were installed at active sites of geomorphic 
change along the muddy coast region in November 1994. The pipes were forced some 3 
to 5 meters into the recent mud deposit, hitting a firm, nearly impenetrable substrate 
(probably the Pleistocene surface, i.e. 6cm to refusal) to assure vertical and lateral 
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stability. Measurements were made from clamps on the pipes to the surface of 
sedimentation. 

2. Sources of Mud and the Deposition Cycle 

The material deposited on the muddy coast originates in the Atchafalaya River 
outflow. This is evidenced by remote sensing imagery (AVHRR, MAS, and MAMS) 
which shows the Atchafalaya River sediment plume extending downdrift from the 
Atchafalaya Bay to the muddy western Louisiana coast. Additional evidence for the 
sediment source is given by the huge quantities of water hyacinth strewn along the muddy 
coast; the Atchafalaya River is virtually the sole source of these plants. Much of the 
subaqueous bottom from the Atchafalaya Bay to the muddy coast is layered with mud 
deposits of varying viscosity. In the regions of greatest shoreface deposition, these 
massive mud deposits have a consistency of heavy cream or “yogurt” and are referred to 
as fluid mud. Sediment cores offshore of the muddy coast typically show several 
centimeters of a light yellowish brown mud overlying a dark gray to black mud of a sub- 
bottom reducing environment. The light yellowish brown sediment is the newly deposited 
fine grained silts and clays that are characteristic of suspended sediment in the Atchafalaya 
River sediment plume. The dark gray to black sediment from beneath the sedimentary 
surface of the inner shelf is episodically exhumed by erosive wave action, creating a black 
“surf zone”. During strong cold fronts, which supply the energy for resuspension and 
deposition events, both fluid mud and the suspended sediment material in the plume are 
sources of deposition material. In the prefrontal phase, onshore winds drive up nearshore 
water levels. If high water levels are attained, the resuspended fluid mud and suspended 
fine grain sediments are washed over the muddy coast and deposited. In the postfrontal 
phase, the cold, dry offshore winds set down the nearshore water levels, stranding and 
dewatering the mud deposits. In the quiescent interludes between strong storms, the 
sediment plume acts to replenish fluid mud by sediments precipitating from the plume. 
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3. Geomorphic Impact Index 

The Geomorphic Impact Index (GI 2 ), which relates atmospheric forcing to coastal 
geomorphology response, was tested using three primary predictors: wind energy, water 
level, and the presence of deposition material (i.e. suspended sediments). The model tests 
the theory that southerly wind energy, in the presence of suspended sediments, translates 
into onshore sediment transport along the Chenier Plain muddy coast of western 
Louisiana. The combination of large wind energy, high water levels, and much available 
fine grained sediment is the most potent condition for coastal geomorphic accretion. This 
set of conditions often occurs in the springtime when southerly winds ahead of cold front 
passages (often stalling over land) appear to be maximized and spring river floods supply 
vast quantities of sediment. Wind energy was monitored using surface observations from 
the Lake Charles, LA NWS station. Optimally, surface winds should be monitored in the 
nearshore waters just offshore; however, the NWS data represented the best continuous 
long term surface wind data set available within close proximity to the coast. Water level 
was represented by the Amerada Hess USACE station in the Atchafalaya Bay, about 80 
km east of the Chenier Plain study site. Sediment availability was monitored using the 
Atchafalaya River stream flow volume data provided by the USACE. It is noteworthy 
that the Atchafalaya River sediment plume/mud stream does not consistently follow the 
prevailing westward drift of the coastal current, but is variously wind driven to the east, 
south or west (as shown by MAMS, MAS, and AVHRR imagery). This phenomenon is 
not accounted for in the use of streamflow data to estimate the presence of depositional 
material along the muddy coast. However, findings of this work indicate that the Chenier 
Plain muddy coast is not sediment starved. 

The investigation relied upon regression techniques. Aerial photography 
observations (2/15/91, 2/28/92, 4/9/93, 2/2/94, 1/8/95, 4/08/96, and 6/3/97) were used to 
manually estimate incremental mud area (mud accumulation/loss on the shoreface from 
year to year) for the muddy coast region of the Chenier Plain. These geomorphic 
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measurements were considered to be accurate to within about 10-20%. Wind energy data 
were summed between the dates of the geomorphic measurements and separated into 
water level bins (-0.5’ to 0.0’, 0.0’ to 0.5’, 0.5’ to 1.0’, 1.5’ to 2.0’, >2.0’). Importantly, 
it was found that wind energy was most correlated to incremental mud area on the Chenier 
Plain when water levels at Amerada Hess were highest (greater than 2 feet). This implies 
that water levels must be near or exceed the elevation of the Chenier Plain muddy 
shoreface to deposit new mud onshore. A separate regression of river discharge and 
incremental mud area showed no identifiable relationship, indicating that fluid mud is 
persistently present along the Chenier Plain (although its replenishment by westward drift 
of the Atchafalaya sediment plume is intermittent). 

The GI 2 expression used is 

GI 2 (n) = -319500 + 28 . 06 £"_ , WE( 2) 

where GI 2 (n) is the Geomorphic Impact Index for date number n. WE(2) is the summed 
Wind Energy when water levels exceed 2 feet. This expression explains 64% of the 
variance in the incremental mud area; the relationship is graphically depicted in Figure 2. 
With the exception of the 1993 sample, the model trends quite well with the observations. 
In spite of reviewing available data, no explanation has yet been found for the 1993 
anomaly. Overall, these results are promising in linking the atmospheric forcing to 
sediment dynamics in the coastal zone; however, they are applicable only to the muddy 
Chenier Plain coast of western Louisiana and are based on a limited, though carefully 
drawn, sample ( 1 data point per year for seven years). 
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Date 

Figure 2. Geomorphic change of the muddy coast portion of the Chenier Plain of western 
Louisiana, observed and predicted by the Geomorphic Impact Index (GI 2 ). Observed 
measurements taken from aerial photography; GI 2 is based on surface wind energy during 
high water level events along the Louisiana coast. The physical interpretation of the GI 2 is 
that wind energy causes resuspension and shoreward transport of suspended sediments 
which, if water levels are high enough, results in shoreface deposition. Incremental Mud 
Area is the change in mud area along the Chenier Plain from one aerial photography date 
to the next. The GI 2 model traces the observations quite well, particularly from 1994 
through 1997. No explanation was found for the anomalous behavior of 1993, a year 
when mud growth was observed, but little wind energy was present to drive suspended 
sediments onshore. 







Figure 3. MAS near IR imagery of Atchafalaya Bay, Louisiana ai full (50 meter) 
resolution (top) and at simulated MODIS 250 meter resolution. Most turbidity features in 
the bay will be detectable with MODIS data. Blurring of narrow inland waterways is 
signilicant. Imagery taken trom April 9 1996 ER-2 overflight of Louisiana coast. 





C. Implications for Applying MODIS data to Coastal Studies 

MODIS will be launched in mid- 1998 on board the AM-1 platform. In its polar 
orbit, MODIS will overfly the Louisiana Gulf coast region twice daily. MODIS contains 
two reflectance channels at 250m spatial resolution, a red channel and a near IR channel. 
These channels will be able to depict fine scale features in turbid water of the Louisiana 
coastal zone (Figure 3) and will be very useful for SSC estimates, coastal circulation 
information, and water types analysis (although the thermal band is available at only 1km 
resolution). After the launch of PM-1 in 2000, it will be possible to identify water motions 
using the 250m bands from AM-1 and PM-1 together. 

D. Remaining Questions 

A number of questions remain concerning the response of sediment distribution, 
transport, and deposition to atmospheric forcing and the resultant impact on the Louisiana 
coastal environment. A few of these are: 

1) Can cold fronts be classified by their associated processes so that geomorphic response 
can be predicted by knowing the cold front type? 

2) What are the details of the geomorphic/sedimentological processes of coastal change 
set in motion by cold fronts and tropical systems? 

3) Do different aspects/phases of the cold front system dictate the response of the various 
types of coastal environments (e g. estuaries, marshes, deltas, etc.)? 

4) What are the nearshore oceanic responses to cold front events (circulation change, 
water level change, wave activity)? Are, individual strong cold fronts more important than 
multiple weak fronts to these responses? 



5) What are the Geomorphic/Sedimentological impacts on the inner shelf (subaqueous 
mud redistribution, replenishing sediment supply, sedimentary structure) by cold fronts? 
Are critical levels of storm strength required before important impact occurs? 

6) What are the marine ecological responses to sediment redistribution, transport and 
deposition (nursery environments, oyster beds, breeding grounds)? 

7) What are the consequences of Earth- Atmosphere System Global Change for the 
creation, maintenance, and vitality of fragile coastal ecosystems (fewer and stronger cold 
fronts, increased tropical activity, rising sea level, increased river discharge)? 
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ABSTRACT 

The Gulf Coast Chenier Plain, a type of coast formed down drift of 
major river deltas, undergoes episodes of erosion and deposition variable in 
space and time. They are controlled by changes in coastal circulation tad 
proximity of delta distributaries. Since 1987, NASA has sponsored several 
studies of coastal geomorphology and sedimentation as drive n by cold front 
passages. A major discovery has been that die Louisiana Chenier Plain is 
rapidly evolving from an erosional coast, threatening coastal properties to a 
depositional one changing into a rapidly prograding, broad, low relief marsh 
plain and mudflat The formerly erosional coast is now being "swamped" with 
fine grained sediment discharged from die nearby, newly invigorated 
Atchafalaya and Wax Lake outlet distributaries of die Mississippi River delta. 

Aircraft and satellite borne sensors have enabled us to follow die progress of 
changes and identify the active processes.. 

L0 Introduction 

Dynamic air-sea-land interactions in micro tidal(l 3-4 S cm) coastal environments 
result in primacy of atmospheric and river discharge driven processes which: (a) power the 
oceanic/estuarine circulation, (b) form various coastal water types, (c) distribute suspended 
and bottom sediments, and (d) drive coastal erosion processes. Sedimentation along the 
Louisiana Chenier Plain coast is in transition from an erosional to a depositional regime. This 
transition has been tri gg ered by die following events: 

1. extension of die active lobe of the Mississippi to the outer shelf edge, resulting in 
decreased flow efficiency where a new channel is now favored, that of die Atchafalaya 
distributary, (see figure 1). 

2. die sedimentary infilling of the Atchafalaya River basin, an 180 km long NNW trending 
topographic low of swamps and lakes filled in with lacustrine deltas(Tye, 1986) and overbank 
river deposits. 

3. The regional coastal circulation with a prevailing westward drift 

The Louisiana-Texas Chenier plain extends from just west of Marsh Isla n d , Louisiana 
(figure 2 ) to Galveston Bay. It is a circa 15km wide coastal plain, consisting of an array 

• Presented at the Second IntemaiionaJ Airborne Remote Sensing Conference end Exhibition. Sea Frmdeoa, California, 24*27 June 1996. 


D-689 





Figure 1. Regional chartlet showing the Mississippi River, die Red River, the Atchafalaya 
River junction at die old river control structure the Atchafalaya Basin and the Mississippi 
River Delta. 


of nearly shore parallel shell/sand beach ridges(Cheniers) and marsh plant stabilized mudflats. 
The beach ridges support long linear stands of large, sturdy live oaks that show up in sharp 
contrast to the marsh grasses. " Chenier" is the french word for "oak", thus the name. The 
Atchafalaya distributary would easily by now have captured the entire flow of the Mississippi 
but the negative economic consequences have prompted the Corps of Engineers to erect a 
control structure allowing 30% of the flow down the Atchafalaya distributary and 70% down 
die present main channel. 
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Figure 2. The Atchafalaya mud stream, it origin, path along the coast and the easternmost 
portion of the Chenier Plain. 


The rapid coastal changes have been "tracked" by a series of: 

1. ER-2 overflights with Modis Airborne SimuIatorfMAS), Multispectral Atmospheric 

Mapping Sensor(MAMS), and color aerial photography. 

2. Leai] et overflights at two altitudes, 6000ft and 39,000ft with the Calibrated Airborne 

Multispectral Scanner(CAMS) plus aerial photography 

3. Surveillance of coastal and continental shelf waters via NOAA satellite AVHRR 

imagery using water temperature and turbidity measurements to track the transport by 
coastal mudstreams. 

4. multiple ground truth field trips. 

This study is focused on a once every 1000 years geological event; die abandonment of one 
major delta lobe(the present birdfoot delta) and die birth of a major new delta lobe by the 
Atchafalaya distributary. . Here, for the first time we observe the earliest phases of a delta 
lobe shift, which involves a massive transfer of water and sediment from one lobe to another 
by "stream capture". It is taking place in "slow motion” due to man's interference, but allows 
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us to follow and document the progress carefully, and to understand these earliest deposits 
These are perhaps most likely to be preserved in the geological record, protected from #.,™ 
by massive later deposits accumulating above and to seaward. S10n 

2.0. THE COASTAL CURRENT AND THE ATCHAFALAYA MUD STREAM 

sunnlv T *T h Coast receives its essential sediment 

supply from the waters and sediments of the Atchafalaya Distributary that are injected into 

wLTaT^A "ST? c CUrTent T1ie rfver out£low « normally a stream of muddy 

^ tchafaI ^ a M n d Stream"(Wells and Kemp, 1981) entrained in the coastal 
current It becomes a discharge "plume" only when it flows offshore during slack current 
periods or when it is driven offshore by Northerly winds. 

2.1. THE ROLE OF COLD FRONT PASSAGES 

ar.™ r® J 0f . ll !f , coas ^ 1 cuma,t “ d coo'. <«*id waters of the Atchafalaya mud 

S 3) C r„M “ AdV “f ed Ve0r Radiometer imagety 

^,Zc= ) of C e: 1 < f,X P S«* re “ nOSP, " riC f ° rCin * P Tocesses that provide for a 

L Pre-frontal strengdiening of southerly winds, increased wave action on the coast 
Mdwntmental shelf, setup of sealevel, spillage of inner shelf waters into estuaries, 
esuspension of sediment and a temporary wanning trend. 

2. Overrunning by a frontal passage squall line, a aune of extreme atmospheric 

turbufcnce, strong and vanahle winds, much oceanic and estuarine vertied mixing, 
sediment resuspension, and rainfall ^ 

T Cold Air outbreak behind the cold front passage. Conditions change rapidly to 
^ong northwesterly wmd^ which are cold and dry, which continue forcing vertical 

sedimenHarf 33 ^ ♦ ^ °r waters > forcing of sealevel downward to spill the 

sediment laden waters out of the estuaries into the inner shelf environment These 

tiTJZt T™ 4e Prevai,in * current from its westerly flow 

% * southeasterly direction. Following most cold fronts a significant SE 

AW?* P ,ume ma V be seen from the mouth of the Atchafalaya Bay 

(figure 3a). (Huh, Rouse, and Walker, 1984,and Huh, et al, 1991) 

a^° timCS per year combined with the high sediment loads of the 

seasonriet^T 11 ^ ° ’^ T **** rangin8 to m >°°° ton * daily in flood 
ason), are coastal processes of fundamental importance to this region. They force 
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Figure 3a. NOAA Satellite Advanced Very High Resolution Radiometer, im age of ihe coastal 
waters. This image is a calculated water leaving reflectance calculated by subtraction 
of AVHRR channel 2 from channel 1 and applying an aerosol correction 
(Stumpf and Pennock, 1989 ). 

3b. NOAA satellite AVHRR thermal image showing the cold coastal current and its 
westward drift 


11-693 







exchanges of water and sediment between the inner shelf and the rivers/estuaries of the 
region. During the sea level setup period they block die river discharge and back these fresh, 
se diment laden waters into the coastal marshes and wedands. When sealevel is set down, 
these waters drain seaward, leaving most of their sediment loads among die roots of the 
marsh grasses which act as filters. This replenishes nutrients of the wetland soils and 
contributes to regional sedimentary upbuilding(Reed, 1989 ). Reed further finds that 
Maximum sedimentation is associated with strong southerly winds both causing increased 
flooding and mobilizing sediment from open bay areas.” 

These processes of river runoff, storm driven exchange, temporary water storage, 
storm surges, important water level rise and fall events, all produce a mosaic of coastal water 
types in this coastal area. Without a reliable salinity tag, we have used water temperature 
and turbidity to quantitatively classify water types and their distribution patterns to identify 
sources(Huh, et al, 1996). We have five recognizable types in the vicinity of the Atchafalaya 
Bay and river mouth: 

1. River discharge waters 

2. continental shelf waters 

3. ambient bay waters 

4. Marsh drainage waters 

5. soil water 

These water types have been detected with die Multi spectral Atmospheric Mapping 
Sensor(MAMS) and die Modis Airborne SimuIator(MAS). They are most clearly 
differentiated just after a cold front passage, readily detectable through die clear diy air of the 
cold-air -outbreak and largely cloudless atmospheric high pressure system. In the quiescent 
period following a cold front passage, the physical differences progressively fade, with solar 
warming of surface temperatures and settling out of suspended sediments. The coastal 
circulation returns to a flow to the west and only die estuarine discharge plume that has been 
driven offshore and the increment of sediment in die marsh remain as evidence of die storm. 

2.2 THE ROLE OF RIVER DISCHARGES 

The Mississippi River is the sixth largest river in the world in discharge with an 
annual average flow of 14,000m j /sec, injecting some 580 km’tyr of fresh water into the Gulf 
of Mexico(Walker, et al., 1994). Springtime river flooding, March-May produces the 
sediment discharge maximum of 300,000 tons daily.. This sediment flood is coincidental 
with the greatest number of cold fronts that stall in die vicinity of die coast, prolonging the 
prefrontal conditions of southeasterly winds, which setup the westward coastal flow, strong 
wave action, and sealevel rise along the coast. The discharged sediments from the 
Atchafalaya river then flow westward toward die Chenier Plain. The fine grained sediment 
concentrations in die river range from 250-400 mg/1, but are found to increase to more than 
800 mg/1 in die mud stream just outside the mouth of Atchafalaya Bay(Wells and Kemp, 
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1981). The concentrations fluctuate with the discharge of the river, which along with the 
storms, modulate the volume of sediment rich waters. The mud stream, bypasses Marsh Island 
and is seen to attach to the eastern-most portion of the Chenier Plain, in die vi cini ty of Fresh 
Water Bayou. Along this coast, the shore face is a mudflat, with a 2 to 20m wide fluid 
mud strip(consistency of heavy cream) just seaward along die coast Offshore of this fluid 
mud zone there is a two layer water column consisting of muddy seawater on top and the 
fluid mud beneath(see figure 5b). The two are distinguishable by the character of the waves 
that propagate through them. The fluid mud based muddy water transforms wind waves to a 
parade of long crested solitary waves, the denser fluid mud zone attenuates the waves 
completely. 

3.0 GEOMORPHIC DEVELOPMENT OF THE PROGRADING CHENIER PLAIN COAST 

Wells and Kemp(1981) noted that "Deposition of fine-grained sediments from the 
Atchafalaya mud stream can be noted by the transitory mudflats that have formed along the 
eastern margin of the chenier plain." This condition continued until 1987, when a dramatic 
changes took place as mud began to accumulate rapidly along this coast Over the years 
between 1987 and 1994, the muddy shore face of the eastern chenier plain began prograding 
rapidly, at rates of approximately 50 m/yr. Low and high altitude aerial photographs show 
the rapid progradation and the geomorphic changes taking place(see figure 4a and b, 5a and b, 
and 6). Figure 4a, representing approximately 1 km of shoreline in winter, 1987 shows the 
wind waves approaching the shell beach face typical of erosional conditions Figure 4b 
shows die same site in 1990, when extensive mud deposits now dominate die shoreface. 

Though averaging approximately 50 m per year progradation into the open Gulf of 
Mexico, this buildup does not occur uniformly along die coast The mud accumulations occur 
in a pattern of two arcuate mud bars, convex seaward along die coast, see figure 7 which 
illustrates one of them. The eastern-most zone of mud a cc umulation is just west of Fresh 
Water Bayou, figure 4. Figures 4, 5, and 6 show the progression of seaward progradation, 
with a transition from a shell beach shoreface(1987) to a zone of mud offshore (1988), to a 
sparsely plant/colonized mudflat(1990), to a half kilometer band of dense, robust, salt tolerant 
marsh (see figure 5b). This growth has been rapid and persistent Two processes are critically 
important to this land building "event", (a) delivery of the fine grained sediments from die 
coastal waters to the shoreface and marshy mudflats directly landward, and (b) stabilization of 
the watery fine grained sediments after deposition. 

3.1 DEPOSITION OF FINE GRAINED SEDIMENT AND GEOMORPHOLOGY OF 
DEPOSITS 

Sediment is delivered to the shore face by sealevel oscillations occurring on a variety 
of scales, from hurricane or winter storm surge, to wave runup. We have observed incident 
solitary waves runup over the nearly horizontal surface of sedimentation at die shore face. 

They runup to a stop, instantaneously drop their sediment load, then yield a clear, sediment 
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Figure 4. (a) coast just west of Freshwater Bayou, about one kilometer stretch. Note well 
developed shell(erosional) beach with windwaves approaching die strand. This was flown in 
Dec., 1987. 

4.(b) Same site in 1990, when extensive development of mudflats at the shoreface. 
marks die transition to an depositional coast due to a flood of fine grained sediment 
transported hear by the coastal current 
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Figure 5. (a) Low altitude(6000ft), 3m resolution CAMS near 1R image of die triple canal 
site 3 km west of Freshwater Bayou. Note the shore-normal in the marsh and the 
linear ponds formed behind die mudbar zone of accretion 


(b) Low altitude oblique of the mudbar, shows from lower left to upper right die 
following features: marsh growth. Shell beach colonized by plants, muddy water of pond 
behind the mudbar, the mudbar with scattered plants and a Ain veneer of fine shell m 
a zone of fluid mud at the shore face(with no waves) and then the offshore two layer water 
column with solitary waves propagating landward. 
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Figure 6. (a) Hie shore face along a 3 km stretch of the coast, with wind waves directly 
impacting die strand, Dec, 1987. 

(b) The newly deposited fluid mud along die coast, January, 1988, no plant 

colonization. r 

(c) The same ate in 1990, with extensive development of the mud flats with linear 

patterns of mud depoation. Plant colonizations on the landward side against the canal endings. 


H-698 









■&!&* 

' "=' .... .= / i .I* . . 

free backwash into the Gulf We have observed deposition on a cm scale, side by side along 
the coast within 10m, with similar scale wave erosion processes. The difference is 
deposition sites are fronted seaward with a local bolus of fluid mud and die eroding area, not 
so protected, is actively eroded by the incident windwaves. After individual cold front 
passages, we have found fresh sheets of gelatinous mud, both widi and without shell sand 
basal layer. In fact the mudflats are built of these multiple sheets of sediment, as can be seen 
when wave erosion of the mudflats occur. Mudbars are formed along die coast leaving linear 
ponds behind the active shore face(see figure Sa and b ). The shore face is convex upward 
in profile when accretion is dominant process and concave upward, when wave erosion is 
taking place. Hurricanes create either an extreme version of the cold front deposit or if the 
predominant winds are off shore, i.e. Hurricane Andrew, they create a widespread sheet of 
fine grained sediments penetrating into die marshy lowlands, with little other discern able 
change. 

3.2 STABILIZATION OF THE DEPOSITS ON THE C HEN I ER PLAIN SURFACE 

The freshly deposited fluid mud is transported landward by the turbulent storm surge 
overriding previous deposits and penetrating in among die roots of the marsh grasses. In die 
backwash, seawater outruns die fluid mud, which is trapped by marsh grasses or stranded on 
die shore face. It has the consistency of "jello" or "tofu* initially, with gravitative water 
drainage out seaward through scattered channels(see figure 7, small shore-normal dialing ^ 
by sheet flow, and by wind driven evaporation. This dewatering is die first phase of the 
consolidation process. The mud cracking process begins cm two scales, by dewatering on 
scale of 20-3 0cm polygons, and then soon by sun creating polygons of a size 

controlled by deposit thickness, generally 8-10 cm on a aide. They progressively d^ccate 
and the mud beach has been converted to a cobble beach wi thin a few days. Wave action, if 
any rolls these cobbles around, rounding them, and when rolled over a shell h«ch become 
what is called die "armored mudballs", sticky cobbles of mud that collect a coating of shell 
fragments. The stabilization process thus involves dewatering, mud cracking, desiccation and 
hardening of mud "cobbles", which now armor the coast; requiring strong wave action to 
mobilize. These changes combined with a zone of fluid mud which strongly attenuates die 
incoming waves clearly demobilizes these deposits. The final stabilization is accomplished by 
die colonization by Spartina and Panicum sp. marsh grasses, which grow robust and well 
developed 2-2.5 m tall, nearly impenetrable on foot with tight interlocking root systems that 
provide strong protection from high wave action. This growth is slack during winter, but 
spring and summer marsh grass growth are vigorous and extensive. Note, along with the 
fluid mud brought by coastal currents is die transport of water hyacinth. These 
plants, coming from the Atchafalaya River, are discharged offshore, die from exposure to salt, 
are carried to the prograding coast by the currents, and are heaved up against die marsh by 
the ton, creating a natural plant fertilization process. The water hyacinth acts as a fer tilizer , 
armors die coast, and further verifies that the sediment source is truly die Atchafalaya River 
and not just onshore drift from die continental shelf. 
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cast over the shore face by storm surges. Length of photo 4 km., 1990. 


4.0 SUMMARY AND CONCLUSIONS 

A major coastal sedimentary and geomoiphic event has been "captured" in a series of 
remote sensing experiments conducted along the Louisiana Coast A transition from an 
erosional to depositional coast is taking place in the eastern portion of the Louisiana chenier 
Plain, still a small portion of the entire Chenier Plain coast This depositional zone is 
propagating westward along the coast at a rate of about a km/yr. Using the AVHRR, 
suspended sediment algorithm, the mud stream of the Atchafalaya River, entrained in the 
westward flowing coastal current has been identified as the source of a flood of fine grained 
sediment in this coastal region. This sediment, deposited on the mudflat shore face by waves 
and storm surges undergo a series of natural consolidation processes which stabilize the 
newly formed land. These deposits have a seaward progradation rate measured to be ~50 
nuyear. The most effective storms in this depositional process are the cold front passages, of 
which some 30-40occur each fall-winter. This is one model of how weather systems of the 
atmosphere, plus coastal circulation of the hydrosphere organize die architecture of 
sedimentary deposits, shaping elements of the lithosphere. Remote sensing with scanning 
radiometers both air and space borne provide essentials on circulation processes, the mud 
stream, and fine structure of die discharge plume that provides for interpretation of processes 
responsible for change. For geomorphology, all possible detail is essential, making high 
resolution aerial photography or very low altitude (< 2000 m) aircraft scanner mission 
essential for definitive measurements. 
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ABSTRACT 

HUH, O.K.; MOELLER, C.C.; MENZEL, W.P.; ROUSE, L.J., Jr., and ROBERTS, H.H., 1997. Remote Sensing of 
Turbid Coastal and Estuarine Waters: A Method of Multispectra] Water-Type Analysis. Journal of Coastal Research , 
12(4), 984-995. Fort Lauderdale (Florida), ISSN 0749-0208. 

Analysis of digital imagery from the 100 m resolution airborne Multispectral Atmospheric Mapping Sensor <MAMS> 
indicates that scatter plots of remotely sensed sea surface temperature versus visible and near infrared subsurface 
reflectance can be used to quantitatively distinguish coastal water types. The Louisiana Gulf coast, a complex region 
of deltas, estuaries, and marshy wetlands, is the setting for this work. Gulf inner shelf and Mississippi River waters 
are the primary source water types, but four additional water types are formed locally and are readily detectable with 
MAMS including: shallow ambient bay water, fresh marsh water drainage, salt marsh drainage, and soil water drain- 
age. In situ measurements of suspended sediment concentration and sea surface temperature in the Atchafalaya and 
adjacent bays have verified the characteristics of these water types. It is found that under the forcing of atmospheric 
cold front passages, water type differentiation is enhanced, creating a mosaic of water types in these coastal waters. 
Multispectral remote sensing of these complex, variably turbid coastal waters, provides data useful for studies of 
water type formation and coastal circulation processes. 

ADDITIONAL INDEX WORDS: Atmospheric cold fronts , suspended sediments, coastal circulation, Multispectral At 
mospheric Mapping System, SPOT images, atmospheric correction. Aircraft multispectral scanner, Chenier Plain, At- 
chafalaya Bay, Atchafalaya River. 


INTRODUCTION 

Mapping the presence, distribution, and boundaries of var- 
ious water types is important to a number of problems com- 
mon to coastal regions around the world. In these environ- 
ments of complex, time dependant air-sea-river-land inter- 
actions, coastal circulation is the essential variable. It is the 
primary control on physical properties, affecting the fate of 
toxic pollutants, the movement of fish larvae, the location of 
fish schools and other marine life, as well as location of ero- 
sion and sedimentation sites. The dynamics of coastal and 
estuarine fronts bounding the water types tend to concen- 
trate water borne pollutants (Klemas, 1980). 

From space and airborne optical imagery, impressive com- 
prehensive views of coastal waters have been obtained for 
many years. Aircraft and manned spacecraft photography, 
and Landsat imagery for example, have provided intriguing 
color imagery of coastal waters. The coastal region often ap- 
pears as a patchwork of plumes, mud streams, and zones of 
various colors and transparencies, distributed from river 
mouths and surf zones to as far away as the continental shelf 
edge. However, for most applications, qualitative depiction of 
color variations is of very limited use. In situ oceanographic 
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analysis traditionally measures temperature and salinity as 
parameters of primary importance. Without remote salinity 
sensing systems at appropriate spatial resolution in the near 
future, water temperature and turbidity /reflectance mea- 
surements provide the useful parameters from airborne and 
spaceborne sensors. 

The purpose of this paper is to demonstrate that moderate 
spatial resolution (100 m or better) multispectral remote 
sensing of surface temperature and turbidity can quantita- 
tively discriminate water types in complex coastal regions 
with high levels of sediment loading. The Multispectral At- 
mospheric Mapping Sensor (MAMS), a multispectral imager 
that flies on board NASA’s ER-2 research aircraft, provides 
visible, near infrared and infrared thermal imagery at 100 m 
resolution. This study is an analysis of MAMS data and how 
it depicts variability of coastal waters. The System Probetoire 
pour Observatione de la Terre (SPOT), a spaceborne multi- 
spectral imager, also provides useful high resolution (20 m) 
visible/near infrared imagery. While the SPOT data are use- 
ful for many inshore studies, the need for a thermal channel 
restricts it’s utility here. 

This study is focused on the Atchafalaya-Chenier Plain sed- 
imentary system (Huh et at ^ 1991) shown in Figure 1. The 
Atchafalaya River brings 152 km 3 of turbid fresh water 
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(Walker and Rouse, 1993) and 88 x 10® tons of sediment 
(Roberts and Van Heerden, 1992) into coastal and Gulf of 
Mexico waters annually (Walker and Rouse, 1993). In this 
deltaic coastal zone, these waters meet in a complex of en- 
vironments including fresh water marshes, salt water marsh- 
es, estuaries, bays, tidal flats, river mouths, deltas and inner 
continental shelf (Figure 1). Suspended sediment concentra- 
tions range up to 1200 mg/1, and more. This coastal region 
has a broad (100-200 km), shallow continental shelf, is 
strongly influenced by river runoff, is microtidal ( < lm, semi- 
diurnal), and is experiencing active subsidence (tectonic-com- 
pactional- fluid withdrawal) of 1.65 cm/yr. (Penland and 
Ramsey, 1990). It will be shown that within this complex of 
environments, atmospherically driven processes alter the 
physical properties and seston content of local marine and 
river waters to generate additional water types, 

MULTISPECTRAL REMOTE SENSING SYSTEMS, 
MAMS AND SPOT 

A series of 2-4 hour NASA ER-2 aircraft missions have 
been flown over coastal Louisiana acquiring MAMS imagery 


and color aerial photography. The MAMS is a cross track, 
line-scanning, multispectral radiometer (Jedlovec et al n 
1989; Moeller et al , 1989). At the operational altitude of 20 
km the 5 milliradian aperture yields 100 m resolution at na- 
dir in 12 visible/infrared channels (Table 1). The swath width 
is 37 km. Calibration of the infrared channels is accomplished 
by viewing two on-board blackbodies of known temperature 
during each scan. Visible and near infrared calibration is 
measured optically in the laboratory before and after MAMS 
deployments through the use of an integrating sphere at 
Ames Research Center, California. Voltage output from each 
MAMS radiation detector is recorded digitally (8 bit preci- 
sion) on board the aircraft and converted to geophysical units 
dynamically in the laboratory. A dynamic calibration soft- 
ware package (Jedlovec et ai, 1989) on the Man-computer- 
Interactive-Data-Access-System (McIDAS, Suomi et a U 
1983) at the University of Wisconsin converts raw digital 
counts to radiances for the visible and near infrared channels 
(channels 2-8) and to radiances and blackbody equivalent 
temperatures for the thermal channels (channels 9-12). 
Earth location and rectification of the imagery is accom- 
plished using McIDAS navigation software in conjunction 


Journal of Coastal Research, Vol. 12, No. 4, 1996 





Huh H al. 


QfiC 

900 


Table 1. Specifications of remote tensing systems 


System 

Spectral Banda 


Spatial Resolution 
Swath Width @ 
Altitude 


C2CS 

Band* 

MAMS 

Band* 

I. 0.43-0.45 *un 

1. 0.42-0.46 pm 

2. 0.51-0.53 

2. 0.45-0.52 

3. 0.54-0.56 

3. 0.52-0.60 

4. 0.57-0.67 

4. 0.86-0.68 

5. 0.60-0.78 

5. 0.70-0.80 

6. 0.66-0.83 

7. 0.72-0.99 

8. 0.83-1.05 

9. 3.47-3.86 
10. 3.47-3.86 

6. 10.5-12.5 

11. 10.55-12.24 

12. 12.32-12.71 

825 m 

100 m 

1,600 km 

37 km 

965 km 

20 km 


Part of 
Spectrum 


1. 0.50-0.59 pjn 

Violet 

Blue 

Green 

2. 0.61-0.68 

Yel/Org 

3. 0.79-0 89 

Org/Red 

Red 

Near rR 

20 m Multupeetrm] 
60 km 

Near IR 
MidIR 
Mid IR 
Therm IR 
Therm IR 


840 km 


with the ER-2 Inertia] Navigation System (INS) data. The 
INS provides aircraft earth position and orientation data re- 
corded at nominal 5 second intervals (Jedlovec et al , 1989) 
For detailed information on the MAMS optical scanning ra- 
diometer see Jedlovec et al. (1989). 

SPOT, the French mapping satellite, carries a “pushbroom 
type imager”, an array of charged couple devices (CCDs) 
which make radiometric measurements in three parts of the 
spectrum, channels HRV-1 (0.50-0.59 pm), 2 (0.61-0.68 pm), 


and 3 (0.79-0.89 pm). It provides multispectral digital im- 
agery at 20 m resolution. In waters of low turbidity HRV-3 
provides excellent land water discrimination. In these turbid 
waters it resolves water types containing concentrations of 
suspended sediment (Figure 2). This system has no thermal 
infrared sensors for remotely sensing sea surface tempera- 
ture. The utility of the MAMS data provides a basis for an- 
ticipating the utility of the MODIS sensors of the Earth Ob- 
serving System when they come on line in the late 90’s. 



upfcrthai.. tart ’SST *» " — « W- wZiU. «, 
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Figure 3. (a) Image location map; (b) MAMS uncorrected band 3 (.56 pm) imagery; (c) after cunglint remora]; (d) after atmospheric scattering removal, 
leaving an image of scaled water leaving reflectance (x * 500- R). Sunglint correction is not attempted over land surfaces. 


ATMOSPHERIC CORRECTION 

In order to isolate the upwelling radiance from within the 
water column, atmospheric contributions to total radiance 
(molecular and aerosol scattering, ozone absorption) in 
MAMS visible and near infrared channels are removed using 
a single scattering model (Moeller et ai„ 1993; Gumley et 
al n 1990). The contribution by specular reflectance (sunglint) 
is estimated using principles outlined in Crackneil (1993), by 
which it is shown that visible sunglint is related to sunglint 
at thermal infrared wavelengths. In this case, 11 pun sunglint 
radiance is assumed small compared to the 3.7 pm sunglint 
radiance. A difference of the 3.7 pm (T >7 ) and IT pm (T u ) 
temperatures is taken as representative of the sunglint at 3.7 
tun and the visible sunglint, for channel k is computed 
by 

Km = a '< T 3-t * T u) + B 

The coefficients A and B are determined for each visible/near 


infrared channel using plots of (T* 7 - T u ) versus visible ra- 
diance in regions of sunglint, and fine-tuned by the user 
through close inspection of the sunglint corrected imagery. 
Effects by atmospheric water vapor absorption and calibra- 
tion uncertainty, as well as dear water radiance are implic- 
itly accounted for in the value of B. Instrument viewing angle 
effects on A and B are not explicitly removed, however, the 
impact appears minimal as no angular biases are evident in 
the MAMS atmospherically and sunglint corrected imagery 
(see example in Figure 3). In this figure, cross track vari- 
ability is evident on the right side of the images 3b and 3c, 
in the relatively uniform inner shelf waters. Figure 3d shows 
a “flat" reflectance field offshore as both sunglint and atmo- 
spheric scattering effects have been removed After correcting 
for sunglint and atmospheric scattering/ absorption effects, 
the resulting water leaving radiance is converted to a sub- 
surface reflectance ratio (water leaving vs incoming radi- 
ance). Hie procedure is outlined in Gumley et aL (1990). The 
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Table 2. Coaetal u**r q,*fy* albt6o m ^g w s (0 65 ^ M ^ g^-NIP) c,. .urfoc, temperature. 


Water Type* 


4 DEC 


5 DEC 


Tamper* 


Temper- 


Gulf water* 

Mean 
St D. 

Bay water* 

Mean 

StD. 

River waters 
Mean 
StD. 

R*®*h water marsh drainage* 
Mean 
StD. 

Salt marsh drain age t 

Mean 

StD. 

Soil Waters 
Mean 
StD. 


aturs 

Albedo* 

•tun 

Albedo* 

ire 

2% 

17*C 

2% 

8.54*C 

0.15*C 

12.6% 

0.1% 

6.8*C 

o.rc 

10% 

0.2% 

11.1*0 

o.orc 

9.3% 

0.2% 

10.1*C 

o&c 

7.6% 

0.3% 

10.2*C 

0.4*C 

8% 

0.3% 

8.5*C 

0TC 

7.2% 

0.2% 

7.9*0 

ore 

6.5% 

0.4% 

6.4’ C 

o.rc 

7.3% 

0.7% 

6.0*C 

1.0*C 

11.2% 

0.8% 

9.rc 

urc 

7.8% 

0.6% 


Description 


tCalculated subsurface reflectance value* (R) converted to albedo by 100*R 
t Differentia ted by location and source 


Warm, low turbidity 
Cold, very trubid 
Cool, moderately turbid 
Cool, very low turbidity 
Cool, very low turbidity 


Vep cold in mid-day, my shallo w wate r environment* with variable 
reflectance, po aai hiy cauaed by aubeqneoua bottom reflectors 


um of a reflectance quantity as opposed to a radiance quan- 
tity is desirable because of varying levels of solar irradiance 
reaching the sea surface, due to changing atmospheric and 
local solar zenith angle conditions. Multiplying the subsur- 
face reflectance value by 100 converts it into a more conven- 
tional albedo value used in Table 2. 

MAMS 11 (im and 12 pm data are used to form a split 
window correction algorithm to remove atmospheric water 
vapor absorption effects on MAMS thermal infrared data 
(Moeller et al, 1993). The resulting sea surface tempera- 
ture (SST) is not corrected for emissivity effects due to sus- 
pended material in the water, the effect of emissivity varia- 
tion on MAMS 11 pm and 12 pm data over turbid waters 
can be reasonably expected to be less than 1 K (Wen-Yao et 
al, 1987). 

METHOD OF MULTISPECTRAL WATER TYPE 
ANALYSIS 

The existence of an inverse correlation between MAMS vis- 
ible'near infrared reflectances and thermal data in coastal 
waters was first pointed out by Moeller et al. (1989). Digital 
interactive inspection of the MAMS subsurface reflectance 
and surface temperature data revealed considerable variabil- 
ity and much useful detail. Greatest discrimination of water 
variability details are possible in MAMS channel 6 (near in- 
frared) and in water surface temperatures calculated from 
Uie MAMS split window channels. Further evidence of the 
importance of the near infrared part of the spectrum for high- 
ly turbid waters can be seen in the comparison of SPOT 


HRV-1 (0.50-0.59 pm vis.) and HRV-3 (0.79-0.89 pm n e a r 
DO data of Atchafalaya Bay (Figures 4a and b). These two 
images are greatly enhanced but not atmospherically cor- 
rected. The faint vertical striping is due to the fact that extra 
enhancement has revealed the very small differences (1-2 
counts) in sensitivity of the Charged Couple Device (CCD) 
detectors of the “push broom" HRV sensor. Figures 5 and 6 
are paired reflectance-temperature MAMS images from 4 and 
5 December, 1990, respectively. The images of both days 
show zones of uniformity with abrupt boundaries between 
a4jacent bodies of water. The flight of 4 December (Figure 5) 
occured in post-frontal conditions approximately 12 hours af- 
ter a cold front passage. Surface winds at New Orleans (MSY) 
Louisiana, were 7.7-10.2 m sec 1 from the north. Water lev- 
els in the Atchafalaya Bay were very low due to the combined 
effects of post-frontal water level set down and low astronom- 
ical tide. The 5 December flight (Figure 6) occurred some 24 
hrs after a slackening of winds. Subsurface reflectance and 
water surface temperatures were extracted from the labeled 
polygons in Figures 5 and 6 and plotted in a tempera- 
turefreflectance scatter plot (Figure 7). This diagram is anal- 
ogous to a temperature/salinity plot used for differentiating 
ocean waters of various origins. The plot quantifies the sep- 
arability of the water types. Knowledge of the local sediments 
and hydrology of the Atchafalaya Bay region along with the 
reflectance and thermal properties in Figures 5 and 6 allow 
us to label each water type in Figure 7. The six water types 
identified are summarized in Table 2. Several of these water 
types (Gulf, Bay, River) have been identified using suspended 
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Figure 6. MAMS images of (a.) band 6 (0.65-0.83 pm) subsurface reflectance and (b.) split window surface temperatures for 5 December 1990, some 36 
hrs after frontal passage. Boxes represent regions of images sampled for 5 December values shown in Figure 7. Note, although several water types are 
apparent on both days, there is a general reduction in contrast of water types from 4 December to 5 Dec em ber, in tandem with a reduction in post-frontal 
winds and a rise in water level (as evidenced by partial submersion of deltaic deposits in the 5 December imagery. 
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^gure 7. Scatter plots of MAMS remotely measured band 6 (.66-. 83 
jun) reflectance and split window SST for 4 December (top) and 5 Decem- 
ber (bottom) 1990. Data sampled from labeled boxea in Figure 5; F » 
Marsh; B * Bay water; S - Soil water, M - Salt Marth; R - 
River. Measurement* of Gulf water* (not plotted) were 17 *C and 0.02 
reapertively. Note the clustering of measurements lor each water type on 
4 December. Some water types also apparent on 5 December, though 
grouping* are more difficult to distinguish. Soil water on 5 December 
lecka true Soil water characteristics but is plotted as Soil water to dem- 
onstrate changing characteristics of colocated waters from 4 D ecember to 
6 December. 


sediment concentration and SST measurements from a lim- 
ited in situ data set collected on Dec 4.J990. It is clear from 
Figure 7 that a bi-spectral approach aids in the separation of 
the water types. For example, the reflectance of bay water 
and soil water are similar, however the temperature signal 
of soU water clearly separates it from bay water. It should be 
pointed out that the “spread” of the soil water reflectances 
may in part be due to bottom reflectance in the very shallow 
water (< 1 meter) near the delta lobes (clear water penetra- 
tion depth for channel 6 of approximately 1 m estimated from 
data provided in Morel and Prieur (1977)) 

Water types plots for other MAMS visible and near infra- 


red channels were produced and examined. Reflectance vari- 
ation among water types was much reduced for the shorter 
wavelength radiation of channels 2 and 3 (blue/green region), 
and for channel 8, near infrared (channel 7 was not used due 
to sensor calibration uncertainty). Also channels 4 and 5 (or- 
angefred) showed better differentiation between reflectances 
of river and fresh marsh waters than channel 6. These spec- 
tral dependencies are subjects of future study. The particular 
characteristics of these water types are dependent on the am- 
bient conditions leading up to and at the time of observation. 
Seasonal variation in suspended material load and sea sur- 
face temperature will shift the plotted points in Figure 6 both 
left-right and up-down. We suspect that the relative position 
of each water type will be generally conserved through the 
cold front season. More observational data are necessary to 
test this idea. 

ATMOSPHERIC FORCING, THE COLD FRONT 
CYCLE 

During the cold front season (October-April), the subtrop- 
ical U.S. Gulf Coast is repeatedly affected by the trailing edg- 
es of passing cold fronts, with mid-latitude cyclones typically 
tracking further north (Palmen and Newton, 1969). Recent 
studies (Roberts et of, 1987; Huh et of , 1991) suggest that 
the annual series of cold front passages have as much or more 
long term impact on these environments than much rarer, 
but stronger hurricanes. Two important facts about cold-front 
p&flsage events are: (1) the frequency of occurrence (30-40 
yr 1 ) and (2) the temporally and spatially ordered pattern of 
change in wind speed, direction, barometric pressure, tem- 
perature and humidity. It is important that the orientation 
of the cold front influences the wind direction, with fronts 
oriented normal (north-south) to the coast having a very dif- 
ferent coastal wind field than those oriented parallel (east- 
west) to the coast (Roberts et of, 1987). 

T^e winter cold front cycle of temperate latitudes is a re- 
peated sequence of events that imparts a burst of kinetic en- 
ergy and an episode of rapid cooling to coastal and near shore 
shelf environments. A detailed study of a strong cold front 
passage and its impact on coastal waters was conducted by 
Huh et aL (1984). From it a simplified descriptive model of 
a cold front passage event is discussed in Roberts et al. 
(1987) and Moeller et al. (1993). A practical consequence of 
a cold front passage is the regional d ominance 0 f a surface 
anticyclone within a few hours or days. The presence of the 
anticyclone optimises optical and infrared remote sensing of 
the surface by typically br in g in g cool, dry, cloud-free air over 
the scene. 

Cold front passages appear to play a significant role in the 
differentiation of water types. The pre-frontal rising sea level 
drives shelf waters over the shore fiscs into coastal marshes 
and impedes river discharge into the Gulf. River water and 
sediment back up into the freshwater marshes and wetlands 
adjacent to the river channels. During the cold air outbreak 
phase, sea level set down induces seaward flow of the fresh 
and salt marsh waters that have been filtered of sediments 
by flowing through the marsh grasses. These low reflectance 
waters are routinely seen off Point Au Fer (salt marsh), Four 
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League Bay, East Cote Blanche Bay, and between the Wpx 
Lake and Atchafalaya Deltas (Figures 1 and 3a). Waterlield 
in the sandy deltaic shoals is also released during sea level 
setdown. The presence of thin streamers of very cold waters 
from the Atchafalaya and Wax Lake Outlet deltas during 
mid-day wanning was the first clue of this soil water drain- 
age. Bay waters, responding to the cold front passage become 
colder and more turbid than the river waters, due to strong 
negative heat fluxes over broad shallow areas and bottom 
sediment resuspension. 

Studies of images and scatter plots indicate that the vari- 
ous coastal waters form, evolve, flow, and generally dissipate 
in time. The few observations we have now through a cold 
front cycle indicate that after the initial burst of post-frontal 
atmospheric forcing, water types are differentiated to a max- 
imum. As the post-frontal forcing weakens (hours to days af- 
ter cold front passage), waters tend to homogenize. This is 
due to processes of cold front forced cooling and settling out 
of suspended sediments. This behavior is apparent in the wa- 
ter types analysis using MAMS data from 5 December 1990 
(36 hours after the cold front passage) in Figure 7. Although 
the 5 December imagery (Figure 6) still shows good variabil- 
ity and indicates the presence of several water types, the plot- 
ted sample sites show how the water masses, especially 
marsh and soil water, have homogenized in both reflectance 
and temperature to the point that it is difficult to separate 
the types from one another. This suggests reduced marsh and 
soil water flow into the bay at the sample sites, a conclusion 
consistent with increased water levels in the bay. (Observa- 
tions at New Orleans showed post-frontal winds had de- 
creased to 2.6-5. 1 m sec 1 from the northeast on 5 December, 
compared to about 7.7-10.2 m sec' 1 from the north on 4 De- 
cember, reducing the water level setdown forcing in the At- 
chafalaya Bay region.) The combination of reduced marsh 
and soil flow with the homogenization of water characteris- 
tics is probably responsible for the lack of separability in the 
5 December plot. Instrumenting key outflow points around 
the Atchafalaya Bay region would further improve under- 
standing of the impact of cold front passages on the cycle of 
marsh inundation and drainage, which is a key mechanism 
by which coastal marshes are supplied with invigorating sed- 
iments and nutrients. 

DISCUSSION 

For coastal oceanographic applications, important ques- 
tions are: (a) what is the physical significance of the temper- 
ature and reflectance measurements acquired by the MAMS 
and (b) how well can this instrument remotely measure 
them? The split window atmospheric correction algorithm al- 
lows a calculated estimate of water surface temperature 
which agrees with in-situ measurements to within 1 °C. How- 
ever these MAMS radiometric thermal measurements origi- 
nate from the top 1 mm of the water column. While this 
leaves the temperature structure of the bulk of the water 
column unsampled, it has been shown that during cold post- 
frontal conditions and under tidal stirring, shallow bay (<1 
to 6 m depth) water temperatures are vertically well mixed 
(Figure 8). Thus, particularly under cold post-frontal condi- 
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Figure 8. A temperature, salinity, and transparency profile obtained 
from a research vessel 24 hours after a cold front passage. The plot shows 
the destratified state of Atchafalaya River Discharge waters in the bay 
in terms of temperature and salinity with a subtle (3%) change in the 
transmittance (suspended sediment concentration) profile. 


tions, one can assume that MAMS SST is representative of 
shallow bay water column temperature. This is not the case 
for suspended sediment load. 

Reflectance of these coastal waters is controlled by seston 
content. In contrast to temperature, the reflectance measure- 
ments originate from varying levels depending on the scat- 
tering/absorption characteristics of the waters as a function 
of wavelength. For highly turbid bay or river discharge water, 
the signal originates from the top few decimeters (Figure 9) 
in most of the MAMS channels. In less turbid waters the 
signal originates from within the top meter (Figure 9). In 
shallow bay environments this can result in bottom reflec- 
tance being detected by the sensor. For this reason, it is ad- 
visable to choose spectral channels in the orange/red to near 
infrared range for reflectance measurements in bay or near- 
shore environments. The ambiguity imparted by detecting 
bottom reflectance possibly contributes to the reduction in 
reflectance range for the different water types observed in 
MAMS channels 2 and 3. In addition, vertical reflectance 
variation in the water columncan he expected (see Nanu and 
Robertson, 1990) due to vertical gradients of suspended 
material. In the case of the Atchafalaya Bay region, vertical 
variation of water types is plausible as marsh, river, and soil 
waters drain into and override ambient bay waters. This also 
contributes to ambiguiety in reflectance observations from 
spectral channels that “see” further into the water column 
(<e.g, MAMS channels 2 and 3). These arguments suggest the 
combined use of orange/red - near infrared spectral data with 
thermal data for identifying water types. 
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(2) the flow of tuspended sediment rich, inshore, marine wa- 
ters through salt marsh areas, with subsequent filtering 
of sediment loed reducing reflectance, 

(3) the resuspension of bottom sediments in shallow waters 
of estuaries and river mouths, through cold front driven 
wave and current action, increasing bay water reflec- 
tance, 

(4) the entrapment of river or estuarine waters as soil water 
within the sand banks of the newly formed deltas. These 
waters drain seaward under very low sea level stands 
associated with the cold air outbreak (post-frontal) phase 
of the cold front cycle. Low temperatures are results of 
the very short term exposure to solar wanning, and 

(5) the cooling or heating of standing, shallow surface waters 
by air-sea beat fluxes, heated or cooled dependant on at- 
mospheric temperatures/humidity, wind speed, water 
depth, and time in the solar diurnal cycle. 


figure 9. Water tranaparency profile* and aeccfai depth measurement* 
in the Atch&f&laya Bay before and after a cold front passage. Theae pro- 
file* show the distribution of light scattering particles (suspended sedi- 
ments) in the water column, in two patterns, (a) pre-frontal pattern, a 
dear layer (75% transmittance) with a bottom turbid layer (50%-0% 
transmittance) of entrained sediment from before the frontal passage, 
and (b) post-frontal pattern, a highly turbid (8%), destratified completely 
mixed layer after the cold front passage. 


TTiese processes, individually or in combination are the pri- 
mary generators of the multiple water types detected in the 
MAMS imageiy. With the aid of cold front cycle forcings, 
these processes are put in motion culminating with an array 
of water types formed or released into the coastal environ- 
ments in the post-frontal phase. 

Study of these coastal waters through analysis of MAMS 
imagery reveals two important general facts: 


Stratified water columns are the norm in the inner shelf 
seaward of the discharging deltas. The active river discharge 
plumes are seen to over-run waters of relict plumes (inactive 
and dissipating) that in turn are overriding warmer, more 
saline shelf waters. These flow westward with the prevailing 
coastal current, a flow interrupted by cold front passages that 
drive episodic eastward flow. The multiple ocean fronts are 
subjects of further studies, 

CONCLUSIONS 

Repeated coverage of coastal Louisiana with airborne and 
spacebome high resolution multispectral scanner data com- 
bined with selective ground truth, is providing a powerful tool 
for coastal research and an important data base for environ- 
mental management. Understanding these complex environ- 
ments has frustrated conventional data acquisition capabili- 
ties for many years. These data are providing an important 
historical record for analysis of long term changes in geo- 
morphology and hydrology of these rapidly evolving Louisi- 
ana coastal-deltaic environments. This record is essential for 
development of environmental management programs impor- 
tant to sustaining and expanding the regional economy. 

Prom the remotely sensed data and surface observations 
we recognize the following coastal water modification pro- 
cesses active in these environments: 

(1) the rerouting of river discharge waters through fresh wa- 
ter coastal marshes by the damming effect of rising sea 
levels associated with the pre-frontal water level setup. 
Filtration of suspended sediment load from river waters 
by passage through grassy marshes renews marsh 
growth and reduces the reflectance of fresh marsh waters 
when they drain in poet-frontal conditions, 


(1) The beet information content on these highly turbid (100- 
1200 mg l' 1 ) coastal waters is contained in the atmos- 
pherically corrected MAMS thermal infrared (sea surface 
temperature from band 11, 10.55-12.24 jim and band 12, 
12.32-12.71 pm) and the orange/red-near infrared 
(MAMS channels 4-6, 0.57-0.67, 0.60-0.73, 0.65-0.83 
jun respectively). The SPOT HRV-3 (0.79-0.89 nm) data 
also gives useful discrimination of water types (Figure 4b) 
at higher spatial resolution. Obtaining time series images 
of these data is important for studies of estuarine and 
coastal pro cesses , 

(2) the complex patterns of temperature and reflectance de- 
fine repeatedly detectable coastal water types, whose 
source, areal extent, and evolution through time can be 
analyzed using sequences of this imageiy. The MAMS 
data contains three basic kinds of information useful in 
identifying water types: temperature, reflectance, and ae- 
rial distribution patterns (i.a, dimensions, orientation, 
boundaries, and associations). Utility of these tools are 
enh an ced when applied in time-series to areally exten- 
sive, complex coastal processes. 

Making meaningful and adequate measurements within the 
regional environment has been too costly or impossible with 
conventional methods. The advent of airborne and space- 
bome sensor systems, global positioning systems, and rapid 
deployment boat operations provide for observations and 
measurements at the necessary time and space scales, and 
under conditions not previously possible. 
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1. INTRODUCTION 

Water motions associated with coastal 
circulation affect many economic, biological, 
ecological, and other concerns. For example, 
timely knowledge of pollutant transport (e.g. oil 
spills) can be critical to minimizing impact on 
endangered marine biology. Along coasts of the 
world, turbid effluent from major river mouths 
contrasts with clear saline nearshore waters. 
GOES-8, providing essentially continuous 
monitoring of North and South America, can 
potentially track nearshore turbidity contrast for 
tracers of water motions. 

From April 1994 into January 1995, 
GOES-8 was positioned at 90°W, providing near 


nadir viewing of turbid waters flowing from the 
Mississippi River distributary system into the Gulf 
of Mexico. Transport of effluent from the 
Atchafalaya River (Figure 1 ). a distributary of the 
Mississippi River, has been under study using the 
MODIS Airborne Simulator (MAS) (King et al, 

1 995) on NASA's ER-2 high altitude aircraft; 
turbid plume response to surface winds has been 
documented (Moeller et al., 1993). In this paper, 
the use of GOES-8 1 km visible data is explored 
as a possfoie source of water motion information 
along the Louisiana coast. GOES-8 1 km and 
MAS 50 meter resolution visible data 
documenting a cold front passage case from 
January 1995 are examined. 



Figure 1 Louisiana coastal region and primary study region for turbid water motions. 
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2. WATER MOTION DETECTION 
APPROACH 

Imagery from GOES-8 can be viewed in 
sequence (looped) using the Man-con^xiter 
Interactive Display System (MclDAS) at the 
University of Wisconsin’s Space Science and 
Engineering Center (SSEC). Researchers at 
SSEC have been active in detecting cloud 
motions in satellite imagery for more than a 
decade (e.g. Merrill et al., 1991). Motions are 
inferred by following trackabie features, or 
targets, from image to image. Targets are tracked 
either manually (human operator visually tracking 
targets with a cursor), or by using a variety of 
automated techniques, most retying on some 
form of cross correlation statistics. For our 
exploratory effort into tracking turbid water 
motions, a manual tracking approach was used. 
Targets are selected at locations where the water 
turbidity gradient is large. 

The accuracy of manually tracked water 
motion vectors is limited by Image to image 
registration acctxacy, image resolution, and the 
operator’s ability to recognize target displacement 
from image to Image. The GOES-8 imager 15 
minute image to image registration performance 
requirement is 1 .5 km and increases to 3 km for a 
90 minute interval (Menzel and Purdom, 1994); 
however, through careful adjustment based on 
identifiable land features, image to image 
registration accuracy can be maintained at a sub- 
pixel level tor small regions (such as the 
Louisiana coast region). For GOES-8 visible data 
sampling resolution of 1 km north-south and 0.8 
km east-west, a sub pixel registration accuracy of 
0.5 km and 0.3 km can be achieved. The relative 
error of the motion vector is the ratio of the 
registration error over the displacement dtotance 
of the vector. Thus, longer time Intervals 
(implying larger displacement) decrease relative 
error, assuming that a given target is accurately 
tracked by the operator over the longer time 
interval. A relative error of 10% requires 5 hour 
and 3 hour time intervals for a water motion 
velocity of 28 cm/s (-1 km/hr) in the north-south 
and east-west directions respectively. Using the 
90 mimAe interval registration performance 
requirement of 3 km, a 5 how time interval results 
in a 60% relative error tor a 28 crrVs north-south 
water motion. Thus, tor typical water motions, 
image registration adjustment reduces the relative 
error of water motion estimates markedly. 


3. RESULTS 

Manual tracking was apptied to hourly 
GOES-8 January 23 Imagery from 1815 UTC 
through 21 1 5 UTC. A cold front had passed 
through the Louisiana coastal region earlier in the 
morning with clearing behind the front reaching 
the Atchafalaya Bay region around 1700 UTC. 
Surface winds of 6 to 8 rrVa from the NNW were 
reported at Lake Charles, LA while 
instrumentation on an ol platform (29° 32.047'N, 
92* 24.814W) just a tow klometers offshore of 
the western Lousiana coast reported winds of 10 
to 12 mto from the NNW. In shatow Louisiana 
coastal waters, post cold front passage winds 
suppress water levels, pushing water essentially 
in the downwind direction (shaltow water Ekman 
flow). Figure 2 shows GOES-8 January 23 vWble 
imagery wfth manualy tracked water motion 
vectors. Relative vector error Is -10-25% for 
these vectors. Whle vectors are sparse, a 
definite indication of water motion in the 
downwind (fraction is present. Water appears to 
be flowing out the mouth of Atchafalaya Bay into 
the Gulf of Mexico. Figure 3 shows the same 
display for January 24, a day In which surface 
winds cHmlnlshed to about 4 mto from the NE (oil 
platform instrumentation) as surface high 
pressure passed to the north of the gulf coast. 
Water motion vectors (15:45-17:45 UTC interval) 
indicate an east to west motion in the nearshore 
coastal waters, essentialy a reversal of the 
January 23 nearshore flow. 

The ER-2 with MAS on board made 3 
overpasses of the Atchafalaya Bay on January 24 
between 15:30 and 17:30 UTC. Water motion 
vectors estimated from MAS imagery are shown 
in Figure 4. Measured MAS registration error 
(through tracking land motions) indicates -10- 
20% relative vector error. It is very evident from 
Figure 4 that many moro targets are identifiable in 
MAS data. These provide a much larger water 
motion vector data base and increase the 
reliability of the estimated water motions through 
opportunky for nearest neighbor comparison. 

MAS shows water flowing from the southeast into 
the western Atchafalaya Bay region, in agreement 
with the few vectors present in the GOES-8 data 
set. The general (fractional pattern shown in the 
MAS data complements that shown in the GOES- 
8 data. GOES-8 vector speeds appear to be 
roughly 25% fester than MAS, a value does to 
the speed uncertainty of both instruments. 




Figure 2. GOES-8 1 km resolution visfcle image of turbkfty pattern along Louisiana coast on January 23, 
1995 (20:18 UTC). GOES-8 water motion vectors (cnVs * 2) are shown pointing in direction towards 
which water is flowing. Water motion towards ESE is apparent on this post-frontal surface wind flow (10- 
12 m/s from NW) day. 



Figure 3. Same as Figure 2 except January 24, 1995 (17:45 UTC). Note reversal of water motions from 
January 23 1 995. Surface winds 2-4 m/s from NE on this day. 



Figure 4. MAS 50 meter resolution viable (0.66pm) imgery on January 24 (1 7:03 UTC). Water motion 
vectors (cm/s * 2) show transport patterns in Atchafalaya Bay region (Figure 1) of Louisiana coast. 


4. SUMMARY 

The contrast of turbid water with dear 
water in nearshore coastal environments for 
producing manual water motion vectors is useful 
with GOES-8 1 km vMbie data. To achieve good 
accuracy (10-20%) for typical water motion 
speeds, it is necessary to Improve GOES-8 pixel 
registration through the use of local landmarks 
present in the imagery. While unable to resolve 
the many small scale (< 1 km) targets present in 
a turbid coastal environment (as shown by 50 
meter resolution MAS data), GOES-8 doee 
demonstrate a capabilty for monitoring general 
water motion patterns, including an apparent 
response to a cold front wind system as shown by 
a reversal of nearshore water motion from 
January 23 to January 24, 1 995. GOES-8 and 
MAS water motion vectors for January 24 
complement each other by displaying transport 
pattern and general speed estimate agreement 
(within the limits of vector accuracy). Manual 
tracking of GOES-8 visible imagery targets wifl be 
investigated in the spring season, when 
distributary outflow is large and nearshore 
turbidity is high. 
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ABSTRACT 

Sediment-rich effluent from the Mississippi River distributary system flows 
into clear Gulf of Mexico waters along the Louisiana coast. The sediment, 
which serves as material for shoreline progradation, is transported by 
nearshore currents to various portions of the coastal zone. Modis Airborne 
Simulator (MAS) 50 meter resolution data from NASA's ER-2 platform is used 
to monitor suspended sediment transport in the Atchafalaya Bay region of the 
Louisiana coastal zone. MAS repeat coverage visual/near infrared images 
from a typical 3-6 hour ER-2 mission, are georeferenced, remapped and 
animated for manual tracking of inwater turbidity features. MAS 
atmospherically corrected red channel (745 nm) data is used to create 
estimates of suspended sediment concentration (SSC) using a regression 
relationship based on in situ SSC samples. The MAS SSC and water motion 
estimates are combined to produce estimates of sediment transport for the 
near-surface water column (upper 1 meter). Results for 24 January 1995, a 
day of light surface winds from the east, show strong transport from the 
Atchafalaya Bay into the East and West Cote Blanche Bays and relatively 
weak transport to the Gulf of Mexico continental shelf. Based on previous 
work, the transport pattern is expected to respond to surface wind forcing, 
especially associated with the cold front wind system. 

1.0 INTRODUCTION 

The Louisiana coastal zone is affected by sediment influx from the Mississippi River 
distributary system. One of those distributaries, the Atchafalaya River, carries high 
sediment loads (up to 1 000 mg/I) into the Atchafalaya Bay region of the coast 
(figure 1), sustaining a nearshore sediment dominated turbid water mass. Sediment is 


* Presented at the Second International Airborne Remote Sensing Conference and 
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Figure 1 . Atchafalaya Bay Region of Louisiana Coastal Zone. 


transported variously into coastal bays, estuaries, and out onto the continental shelf of the 
Gulf of Mexico by variable coastal currents, which are modulated by surface winds and 
astronomical tides. Previous work has identified new (beginning In 1987) sediment 
accumulation on coastal shorelines downdrift of the Atchafalaya Bay region (Huh et al., 
1991). Moeller et al. (1993) and earlier, Roberts et al. (1987) have indicated that the 
pattern of sediment transport is strongly affected by surface winds associated with cold 
front systems. These winds act to drive up and set down water levels in the microtidal 
Louisiana coastal zone. This paper will focus on the transport patterns on 24 January 
1995 as seen by the MODIS Airborne Simulator (MAS) (king et al., 1996). MAS flies on 
NASA’s ER-2 high altitude aircraft, collecting 50 meter resolution data across its 37 km 
scanning swath in 25 laboratory calibrated visible/near infrared and 25 onboard calibrated 
infrared channels (Table 1). MAS and , earlier, the MAMS (Multispectral Atmospheric 
Mapping Sensor) instrument have flown some 20 missions over the Louisiana coast since 
1 988. Three MAS overpasses of the Atchafalaya Bay over a 2 hour period are used to 
generate water motion vectors using georeferenced visible/near infrared data; MAS 
atmospherically corrected visible/near infrared data are investigated to estimate suspended 
sediment concentration (SSC). These quantities are combined to produce estimates of 
sediment transport. 


2.0 DATA SET COLLECTION 

Under dear skies on 8 and 24 January 1995, MAS was flown on the ER-2 over the 
Louisiana coast. On both days, boat teams from the Coastal Studies Institute (CSI), 
Louisiana State University (LSU) collected In situ water samples (from top 1 meter of the 
water column) in the Fourieague Bay region and off the Chenier Plain coast in 




Table I . Spectral characteristics of the 50-chanrtel MODIS Airborne Simulator (MAS), 



coordination with ER-2 overpasses. Locations of in situ data collection are chosen to 
maximize turbidity variance, improving the value of the limited sample size (- 30 samples). 
On 24 January, the ER-2 flew a racetrack pattern, making 3 overpasses of the Atchafalaya 
Bay between about 1 525 and 1715 UTC (45 minute repeat cycle) (figure 2). Morning data 
collection was chosen to minimize sunglint in the imagery. Surface winds were relatively 
light (5 m/s ENE) as atmospheric high pressure lay just to the north. The dear, dry, post- 
frontal atmospheric conditions optimized the MAS data depiction of coastal water features. 



Figure 2. GOES-8 Image on 24 January 1995 with ER-2 Flight Pattern. 



























3.0 DATA PRODUCTS 


3.1 Water Motion Vectors From MAS Data 

Tracking water (and atmospheric) features for estimating motion vectors is not new 

,n t Breakar at al * 1994 )- Historically, both manual (analyst in the production 
stream) and automated (e.g. maximum cross correlation) techniques have been applied. 

f < ?L!!ri in ? at, r 9 water motions wWh MAS - a simple manual approach was applied to track 
turbidity features in MAS orange-red (657 nm) and near infrared (869 nm) data. In a 
manual approach, the analyst subjectively identifies trackable features, determining the 
feature position in each subsequent image. The MAS imagery are georeferenced and 

CO T mon P ro i ection t0 facilitate viewing a temporal sequence of images with 
mal land motion or jitter' (during ER-2 flight, data on the three axis (pitch, yaw, roll) 
aircraft orientation, altitude and nadir latitude/longitude is recorded for use in positioning 

0r l the ea ™ Mewing three or more images in sequence is desirable 
as the addition^ image(s) help to fix the feature’s progress in the mind of the analyst as 
wen as establish the conservation characteristic of the feature. For the 24 January data set 
in tnis paper a three image sequence was animated on the University of Wisconsin’s Man- 

fh^mo Gr [i > e ^ a ? t ve , Data Access System (MclDAS). Velocity vectors are computed from 
the manually determined position of the feature and associated time data. In the future 
automated feature tracking techniques, including feature rotation and deformation (e q ’ 
Kuo and Yan, 1994) will be tested. 

3.2 Suspended Sediment Estimate 

Isolating the inwater scattered radiance contribution requires the removal of 
atmospheric scattering (molecular and aerosol) and sunglint contributions to the sensor 
observed total radiance. MAS visible/near infrared calibrated radiances are 
atmospherically corrected using the approach outlined in Moeller et al. (1993) This 
abroach « based on single scattering atmospheric modeling work developed by Gordon 
i!mll^ fU ^ r ! d i ,ymany(e - 9 ^ ^°rdonetal.. 1983, Guzzi etal., 1987). Removing the 

** water teavin 9 Glance (radiance upwelfing out of the 
water column) quantity which is convertible to a SubSurface Reflectance (SSR) quantity 

< 199 p)- Th® SSR is essentially the water leaving radiance normalized by 
me solar incoming radiance entering the water column. The 24 January MAS SSR for 
channels 1 -9 (547 - 947 nm) were regressed against the fifteen in situ SSC samples 
collected in Rourieague Bay. The best relationships were produced by red to near infrared 
(MAS f* nnels 4 * €: 745 ’ 786 ‘ 827 nm >- This Is not surprising because the water 
!H5 # r ^ fllon of *" s * )e< * rum Is small (- 1 meter or tan in turbid water), 
S* as such the turbidity signal will originate from a level similar to the in situ samples. 

S oiler wavelength channels integrate over larger columns and may pick up bottom 

m , me Sf ! a,low tey waters ' while k>n 9 er wavelength channels see progressively 
l«ss of the water column and are adversely affected by atmospheric water vapor 
absorption in the 900 - 1000 nm spectral region. The channel 4 (745 nm) regression was 

(fi ? UrB 3) * A quadratic provided a good fit of the data, with an RMS 
of 37.8 mg/I. The 2nd order nature of me SSR-SSC relationship is suggestive of the 



i -» 





Figure 3. MAS Atmospherically Corrected 745 nm Relationship With In Situ SSC. 


nonlinear relationship that exists between optical depth and reflectance in the radiative 
transfer equation. 

3.3 Sediment Transport Estimate 

The MAS water motion vectors and SSC estimates are combined with water 
penetration depth to produce estimates of sediment transport ST, 

ST(kg/s) = V * SSC * 2 * D ( 1 ) 

where V is the water motion velocity (m/s), SSC is the MAS regressed suspended sediment 
concentration (mg/I), Z is the water penetration depth (m) and D is the MAS footprint width 
(m). Physically, eqn. (1 ) represents a sediment transport through the cross section 
represented by the MAS footprint width and the water column layer from which the 
sediment scattering signal originates. From work performed by Gordon (1975), Z Is 
estimated to be about 1 meter or less for the MAS red channel at 745 nm. Much of the 




Atchafalaya Bay is about a meter or so in depth; surrounding bays (East, West Cote 
Blanehe, etc.) are probably somewhat deeper although good bathymetry for the region is 
difficult to obtain given the dynamic impact that the sediment influx and 
resuspensiorVprecipitation processes have on the local bathymetry. As such, the sediment 
transport estimates reported herein do not necessarily represent full water column 
transport. However, it is plausible that water motions are fairly uniform in shallow columns 
(with exception of surface drag at subaqueous bottom), and transmissometer profiles in the 
Atchafalaya Bay region can be used to model the vertical profile of SSC (they indicate an 
SSC maximum near the subaqueous bottom which is likely a result of resuspension 
activity). These assumptions, coupled with good knowledge of the water column depth 
could be applied to adjust the MAS sediment transport estimates to full column estimates. 

4.0 RESULTS AND DISCUSSION 

Figures 4 and 5 show water motion velocity and sediment transport data for 24 
January 1995, respectively, plotted over MAS remapped SSC imagery. Vector scaling is 
shown in the figures; vectors point in the direction of flow. For sediment transport 
estimation, the SSC quantity is averaged over a 3 X 3 pixel box centered at the water 
motion vector location (base of vector). Water motion velocity relative error (georeference 
error over the displacement distance) was estimated in the 10-20% range based on 
tracking land motion in the MAS images. This error translates directly into a 10-20% 
relative error in the sediment transport estimates as well. It is dear when comparing 
figures 4 and 5 that relatively high water motion velocity does not imply high sediment 
transport. It’s noteworthy though that the SSC maximum (southwest side of the 
Atchafalaya Bay) exists in a region of high water motion velocity. The high velodtfes may 
be enhandng resuspension activity near this SSC maximum. Outdated bathymetry charts 
indicate that water depths are about 1 -2 meters around the SSC maximum and increase 
towards the southwest (where velocities remain high but SSC is much lower). A review of 
MAS imagery collected on 8 January 1 995 did not show an SSC maximum at this location* 
thus it is not likely that this maximum is caused by bottom reflectance. Continuous 
sediment influx into the Atchafalaya Bay takes place primarily at two locations, the 
Atchafalaya River and Wax Lake Outlet (located at the right and left newly forming deltas, 
respectively). A 26.7 kg/s sediment transport is estimated at the mouth of the Atchafalaya 
River (no trackable feature was available for estimate at mouth of Wax Lake Outlet). Using 
eqn. (1) and assuming constant water motion velocity and SSC with depth, total sediment 
transport from the mouth of the Atchafalaya River (-1200 meters wide, dredged to about 
15 meters deep for shipping purposes) is estimated at 9612 kg/s. This rate comes to about 
3.0e1 1 kg/yr. Approximating the volume of the high turtoldlty zone in the Atchafalaya Bay 
at 20 km X 20 km X 1 meter depth, and using an average MAS SSC of 200 mg/I over this 
domain, the total sediment load of the Atchafalaya Bay Is estimated at 80,000,000 kg. 

Thus, for the estimated total transport rate at the mouth of the Atchafalaya River, the time 
required for the Atchafalaya River to replace the sediment load is on the order of a few 
hours. This simple approximation shows that the sediment transport rate in the 
Atchafalaya River has a continuous significant impact on the sediment load in the 
Atchafalaya Bay. Reducing or increasing the sediment transport in the Atchafalaya River 




Figure 5. Sediment Transport in the Atchafalaya Bay Region on 24 January 1995. 





(e.g. local heavy rains, upstream heavy rain or seasonal drought/flooding) will cause near 
term variation in the Atchafalaya Bay sediment load (tempered by 
resuspension/precipitation activity). This sediment load is available for transport to other 
bays and shoreface. For example, in figure 5, East Cote Blanche Bay Is the recipient of 
strong sediment transport from the Atchafalaya Bay as water velocities in excess of 30 
cm/s advect highly turbid water. Water motion velocities and sediment transport in East 
Cote Blanche Bay are smaller, suggesting a net positive flux of sediment in that region. 

The Fourleague Bay basin also appears to be experiencing a net po sit ive fl ux of sediment 
from the Atchafalaya Bay. No sediment transport from Fourleague Bay into the Gulf of 
Mexico is indicated through the thin southern outlet connecting them. Instead sediment 
appears to be transported into waterways on the east side of Fourleague Bay. 

Previous MAS (and MAMS) data have indicated that local transport patterns vary 
greatly in the Atchafalaya Bay region. This has been inferred from variant positioning of 
SSC maxima and minima as well as through animation of the imagery. SSC maxima have 
been variously positioned Gulfward over the continental shelf, westward in East and West 
Cote Blanche Bay and Vermilion Bay, as well as eastward in Fourleague Bay. Obvious 
indications of SSC transport from Fourleague Bay into the Gulf of Mexico have also been 
observed. SSC patterns are related to variation of sediment transport from the Atchafalaya 
River, resuspension activity, and nearshore transport currents (a function of astronomical 
tides, and atmospheric forcing). A case study in Moeller et al. (1993) showed that 
southwesterly surface winds retarded sediment transport from the Atchafalaya Bay 
westward to the Chenier Plain region. AVHRR data have shown sediment transported 
southward from the Atchafalaya Bay region onto the Gulf of Mexico continental shelf 
following cold front passages (Walker et al., 1992). Most recently, Moeller et al. (1996) 
used GOES-8 1 km visible data to show variation in general transport patterns associated 
with changing surface wind conditions. These investigations point to surface wind as an 
important influence on transport patterns in and around the Atchafalaya Bay region. In the 
future, additional MAS data sets will be collected and analyzed for transport patterns and 
their relationship to atmospheric forcing. In situ data will be expanded to indude 
continuous monitoring of inwater transmission (for SSC data collection enhancement) and 
the deployment of drifting buoys for estimating water motion velocities. 

5.0 SUMMARY 

A demonstration using ER-2 mounted MAS 50 meter resolution data for estimating 
water motion velocity and sediment transport in turbid coastal waters has been presented. 
MAS geo referenced visible/near infrared data viewed in sequence is used to manually 
track turbidity features to yield small scale water motion velocity estimates. The ER-2 
platform shows itself to be of sufficient stability during repeat flight track coverage to allow 
the production of useful water motion velocities from MAS data (within 10-20% relative 
accuracy). MAS atmospherically Corrected visible/near infrared data have been used to 
estimate suspended sediment concentration (SSC) using a regression relationship 
generated with in situ SSC data collected during ER-2 overflight. The MAS red channel 
(745 nm) was chosen as producing the best relationship. Sediment transport estimates 
generated from the wate r motion velocities and SSC data indicate variable sediment 



accumulation in the Atchafalaya Bay coastal region. Such data are informational for 
hydrodynamic, coastal geomorphic, and marine ecological investigations. The relationship 
of sediment transport to surface wind forcing, especially associated with cold front wind 
systems is being solidified by this work and will continue to be pursued. 
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1.0 INTRODUCTION 

The MODerate resolution Imaging Spectrometer (MODIS) is scheduled for launch into 
polar orbit in summer 1998. MODIS has 36 channels in the visible and infrared region at various 
spatial resolutions of 250, 500, and 1000 meters. MODIS, as part of the Earth Observing System 
(EOS), will provide data for generation of atmospheric (e.g. cloud properties, aerosols), oceanic 
(e.g. ocean productivity, sea surface temperature) and land surface (e.g. land surface temperature, 
vegetation index) products. Applications of MODIS data to coastal studies merit further study. 
The MODIS 250 meter resolution visible channels will provide a high resolution map of coastal 
environments on a repeat cycle similar to current AVHRR monitoring. 


TABLE 1. Spectral characteristics of the 50-channel MODIS Airborne Simulator (MAS), 
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The MODIS Airborne Simulator (MAS) (King et al.,1996) was developed as an airborne 
simulator for the MODIS instrument. MAS contains visible and infrared bands that will be 
available from MODIS (Table 1). On NASA's ER-2 platform (nominal 20 km altitude), MAS 
collects 50 meter resolution data across a 37 km swath. Beginning in 1994, MAS has flown 
missions over Louisiana's Gulf coast, mapping suspended sediment transport, nearshore water 
temperature, coastal water types, and circulation patterns in response to cold front wind system 
forcing (Moeller et al.,1993). MAS (and MODIS) applications to coastal harards (e.g. oil spill, 
storm impact) are the topic of this discussion. 

2.0 MAS APPLICATIONS 

2.1 MAS WATER MOTION DETECTION 

On the aircraft platform, MAS is capable of tracking inwater features (e.g. turbidity, 
pollutants, watercraft) in coastal waters through repeated mapping of flight lines over short time 
intervals (hours). This approach has been used for tracking sediment transport in nearshore 
coastal waters (Mueller et al.,1996), an example of which is shown in Figure 1. In Figure 1, 



Figure 1 . Water Motion Velocity Estimates from MAS 24 January 1995 Data. 
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turbidity contrast made it possible to subjectively select and track features for water motion 
estimation. Accurately georeferendng the MAS data is critical, as displacements of recognizable 
wato- features are typically small over the three to six hours that the ER-2 aircraft can spend 
monitoring a given location (note: the development of Unmanned Aerospace Vehicles (UAV) ca n 
significantly extend onsite duration). The presence of land features in the data of Figure 1 
facilitates good georeferendng of the data; GPS systems on aircraft hdp for water-only scenes. 
MAS water motion vectors provide a picture of the circulation responsible for near-surface 
transport of inwater material. 

One possible application, directly tracking oil slick features, has not been attempted with 
MAS data although oil slicks have been detected in visible spectral data when sunglint is present. 
This approach limits direct oil slick detection to daytime conditions. Spectral si gnatur es of oil 
(and other pollutants) need to be identified for nighttime detection. A combined MAS and 
Synthetic Aperture Radar (SAR) data set should be collected to explore syner gism between these 
data types for surfactant detection. Knowledge gained from such a data set would likely be 
directly applicable to MODIS data. 

2.2 MAS WATER TYPE DISCRIMINATION 

Analyses of MAS visible, near infrared (NIR) and longwave infrared (LWIR) data have 
shown that the various water sources, or types, present in the coastal environment can be 
discriminated from one another using multispectral analysis. Diagrams of turbidity versus water 
temperature separate water types. This tool has been demonstrated using MAS data along the 
Louisiana coast (Huh et al., 1996). Turbidity estimates are taken from MAS visible/NIR 
atmospherically corrected data; sea surface temperature estimates are made using a split window 
technique with MAS 1 1 and 12 micron data. Investigations have shown that a myriad of water 
types (saline Gulf, river, fresh marsh, soil drainage from sub-aerially exposed deltaic sand banks, 
bay, and salt marsh) are present after a cold front passage, due to local surfece winds driving 
down water levels in the coastal environment. The presence of these water types is a si gnal of 
ongoing vitality of these specialized environments, especially in the wake of damagin g tropi ca l 
systems. 


Land surface types can also be discriminated using multispectral analysis, ft is anticipated 
that marsh, bare soil/sand, vegetated, burned, etc. regions are identifiable using MAS data Such 
analyses bring additional information to storm impact assessment, including overwash, marsh 
destruction, vegetation and agricultural destruction. Large (hundreds of meters) changes in 
coastal landform (erosion and deposition) are also identifiable with MAS d?ta 

2.3 COASTAL LANDFORM IMPACT 

The Louisiana coast is experiencing gradual but continuous subsidence (from 1.5 to 5 cm/ 
yr) of its landform. Without continued influx of sediment from the Mississippi River distributary 
system, coastal landform eventually subsides below sea level and is flooded. The end result is salt 
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water encroachment and destruction of coastal freshwater environments. Petroleum industry 
interests in the coastal zone require access channeling which further reduces sediment dispersal 
into freshwater marsh environments. Recently, a plan has been enacted to improve sediment 
dispersal along the coast by opening crevasses (breeches in the levee). Monitoring sediment 
dispersal around crevasses is valuable for understanding the efficacy of this manag ement 
approach. MAS has the high spatial resolution and multispectral capability to contribute to the 
assessment. The impact of tropical storms on crevasses must also be monitored. The disab l i ng or 
destruction of crevasses limits sediment dispersal and effectively accelerates landloss. MAS water 
motion estimation provides a dynamic assessment of sediment dispersal in and around crevasses. 

3.0 DISCUSSION AND SUMMARY 

Coastal environments are among the most sensitive of the Earth system. Because of their 
ability to support a wide variety of marine and terrestrial life, they are also among the most vital. 
Management strategies which support fragile coastal ecosystems and yet maintain fisheries, 
petroleum, recreational, and human habitation interests require detailed knowledge of natural and 
anthropogenic impacts on these environments. These impacts include surfactant discharge, 
tropical and extra-tropical storm impact, and man-induced landform change, among others. 

As part of an ongoing effort, data from the MODIS Airborne Simulator (MAS) has been 
used to monitor the Louisiana coastal environment, primarily for suspended sediment dispersal 
and water type identification. From NASA's ER-2 platform, MAS 50m resolution data can be 
collected over 3 to 6 hour periods for monitoring coastal environments with visible, near infrared 
and ther mal infrared spectral data. Applications developed for MAS data will be transferred to 
NASA's MODerate resolution Imaging Spectrometer (MODIS) scheduled for l au n c h into polar 
orbit in summer 1998 (AM platform). MODIS, as part of the Earth Observing System (EOS), will 
play a prominent role in the Earth system science of the Mission To Planet Earth (MTPE) 
program. 

MAS should be viewed as a tool for monitoring water motions and transport during 
events which require knowledge of transport for planning remedial action (such as oil spills). 

MAS is also a tool for assessing impact on coastal environments by storm systems through 
identifying water types and landform change. As a visible/infrared spectrometer, MAS 
measurements are affected by clouds. Coupling MAS with SAR imagery should provide an 
enhanced description of coastal zone conditions. Use of SAR technology to characterize wave 
conditions during storms may provide an important link between atmospheric forcing and sea 
state response. 
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